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Abstract 

We describe the Hamiltonian structures, including the Poisson 
brackets and Hamiltonians, for free boundary problems for incom- 
pressible fluid flows with vorticity. The Hamiltonian structure is used 
to obtain variational principles for stationary gravity waves both for 
irrotational flows as well as flows with vorticity. 



1 Introduction 

In 1933 Friedrichs [9] proposed the functional 

= JJJ [{V^l;Y + v'ix,y)] d 



where ip is the stream function for an incompressible flow, as a variational 
method of obtaining solutions to free boundary value problems. Critical 
points of J are harmonic functions which satisfy the condition 

on the free boundary, given hy if) = 1. The free boundary condition relevant 
to theory of gravity waves, however, is the Bernoulli equation 

V gQ = constant, 



where if is either the velocity potential for irrotational flow, or the stream 
function in the case of flows with vorticity. Thus some other variational 
principle is needed for the study of gravity waves. 

Recently, a variational principle for gravity waves with vorticity was given 
by Constantin et. al. [7], using a direct, "hands on" approach. More gener- 
ally, a variational principle for a stationary wave may be obtained for systems 
possessing a Hamiltonian structure by minimizing the Hamiltonian computed 
in a Galilean frame moving with the wave. We illustrate that approach in 
this study. 

We begin with a brief review of Euler's equations of incompressible flows, 
and the associated free boundary value problems; in §3 we describe the Hamil- 
tonian structure of these problems, for irrotational flows and flows with vor- 
ticity, as given by Lewis et.al. [15]. All the functions under consideration in 
this article, including the free boundaries, are assumed to be smooth. 



2 Incompressible fluid flows 

Let the velocity field of an incompressible fluid in a fixed region D be denoted 
by V. The incompressibility of the fluid is expressed by the condition div v = 
0. We must have v ■ i/ = on the boundary of T>, where v is the outward 
unit normal at the boundary. Euler's equations of motion for the flow of an 
inviscid, incompressible fluid are 

^ = V, = p(vi + (v ■ V)v) = - V(p - g ■ x), (2.1) 

div V = 0, 

where p is the density, p is the hydrodynamic pressure, and g ■ x is the 
gravitational potential. Henceforth we take p = 1. 

Given a manifold P G with smooth boundary, we denote by C^{T>) 
the Hilbert space of vector fields on V with the inner product 

(v, w) = jjj v-w(i^x. 

We denote by Ct, the closed subspace of C'^iV) generated by vector fields; 
of the form w = V p for some function p with finite Dirichlet norm. The 
orthogonal complement = is the space of all vector fields v for which 
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(v, V p) = for all p G W^''^{V). By applying the Gauss divergence theorem, 
we see that if v G C^(T>) and is smooth, say C^, then div v = and v ■ z/ = 
on S = dV, where u denotes the outward unit normal on S. The Hilbert 
space Ca is the space of weakly divergence-free vector fields. We denote the 
orthogonal projections onto £o- and by and Pt, respectively. 

In many applications the fluid is not confined to a fixed region, but instead 
carries the region with it. In such cases, the region T> occupied by the fluid 
must also be determined. Such problems are called free boundary problems 
and occupy a substantial part of the literature on incompressible flows. 

Given an irrotational flow (curl v = 0) on a simply connected domain, 
there is velocity potential (f for which v = V The velocity potential is 
defined only up to an arbitrary function of time; the transformation ip i-^ 
ip + k{t) is called a gauge transformation, and will play a role in what follows. 

The equation div v = implies that if is harmonic. Substituting v = V 
into the second equation in (2.1) we obtain 



for some function of time, which can be eliminated by a gauge transformation 
of the velocity potential. We always choose the gauge to be such that 



everywhere in the fluid. 

An interface between the fluid and another medium, for example air, is 
called a free surface. If the pressure is constant in the air, then it is also 
constant at the surface of the fluid, and we may normalize the pressure to 
be zero at the free surface. Hence we obtain Bernoulli's equation 



where gz is the gravitational potential on the free surface. 

The free surface is given in space-time by </> = 0, where (j){x,y,z,t) = 
z — ({x,y,t). The free surface moves with the fluid, hence the material 




hence 



^t + -{'^vf + gz + p = k{t) 



^t + -(yv? + gz + P = o 



^t + -{'^vY + gz = o, 
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derivative of vanishes, and 
hence 

This is called the kinematic condition on the free surface. 

This collection of equations for gravity waves on a free surface is known as 
Euler's equations for waves on the surface of an inviscid, incompressible fluid 
with irrotational flow in the region V = {{x,y, z) : < z < h + ({x, y, t)}. 
They are 

Aif = 0<y<h + C, 

Ct + ^xCx + 'fyCy = on S] (2.2) 

<^t + ^|V<^P + (7Z = on^; (2.3) 



ifz = on z =0. 

Here, if is the velocity potential of the flow, and ({x, y, t) the displacement 
of the fluid surface from equilibrium. We have neglected surface tension. The 
second equation is known as the kinematic equation; the third equation is 
Bernoulli's equation. At rest, the fluid lies in the region < z < h; g is the 
acceleration due to gravity. The free surface is denoted by S* = {{x,y,z) : 
z = h + C{x,y,t)}. 

The two physical constants in the theory are g and h. Let c denote a 
characteristic velocity (e.g. the velocity of a gravity wave); then /i/c is a 
characteristic time. We introduce dimensionless variables 

{x,y, z) = h{x' ,y' , z'), t = ht'/c, if = chip' . 

The equation (2.3) now becomes 

^;' + ^(VVr + AC = 0, A = ^, (2.4) 

where A is the inverse square of the Froude number. The other equations in 
Euler's system are unchanged under the rescaling. From now on we drop the 
primes and understand that we are working in non-dimensional variables. 
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The Euler equations are invariant under the one-parameter subgroup of 
Gahlean boosts along the x axis, given by 

{x', y', t') = (x - ct, y, t), v'(x', y', t') = v(x, y, t) - (c, 0). 

The velocity potential, however, is determined only up to a function of time. 
Thus the Galilean boosts on the velocity potential are given by 

y', t') = ip{x, y, t)-cx + q{t). 

Under these Galilean boosts, 

^+2(V^) =^ + 2^^^) +^W-2^- (2.5) 
The result follows by direct calculation, noting that 

d d d d d 

dt' dt dx' dx' dx 

Proposition 2.1. Suppose the solutions of Euler 's equations are stationary 
in a Galilean frame moving with speed c. Then Q = '^'f = 0; and, choosing 
q{t) = c^t, the conditions on the free surface are (dropping the primes) 

Proof. By (2.5) the Bernoulli equation in the moving frame is 

As X ^ ±oo ^ while (V'v?')^ ~^ ■ Moreover, y?^, = by the assump- 
tion of stationarity. These conditions force the choice g' = c^, and the result 
follows. The kinematic equation in the moving frame is immediate. □ 

Proposition 2.2. Let v be a divergence-free vector field in a domain T). 
There is a unique orthogonal decomposition, known as the Weyl-Hodge de- 
composition, 

v = w + V(^, (2.7) 

A(p = 0, = V ■ z/; div w = 0, w ■ = 0. (2.8) 
The proof is left to the reader. 
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3 Poisson structures 



Let M be a C°° manifold of dimension n, and let F,G & C°°{M). A bilinear 
form {F, G} is said to be a Poisson bracket if 

. {F,G} = -{G,F}; 

. {{F, G}, i/} + {{G, ff}, F} + {{ff, F}, G} = 0; 

The second property implies that the Poisson bracket is a derivation in each 
of its entries. Hence any H G G°°(M) generates a vector field Xh-, called a 
Hamiltonian vector field on M, defined by X^F = {H, F}. The Hamiltonian 
vector field Xh generates a flow on M; if are a set of local coordinates on 
M, then the time evolution of the on that chart is given by the ordinary 
differential equations 

x' = {H,x'}. 

Due to the fact that the bracket acts as a derivation on each of its entries, 
we may represent a Poisson bracket in the form 

j,k=i 

where W^^{x) is a skew-symmetric matrix. 

If det W ^ then it is easily seen that n must be even. A classical 
theorem of Darboux states that in this case it is always possible to find a 
set of local coordinates, called canonical coordinates q^,p\ (1 < i < ^f^) 
which the Poisson brackets take the form 



" dF dC OFdG 



A manifold with a Poisson bracket is called a Poisson manifold; if the 
brackets are non-degenerate, the manifold is called a symplectic manifold. 
On a symplectic manifold, the Hamiltonian flow takes the form 

■ i_dH_ 

^ dp"^ ' ^ dq'^ 
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In this paper we shall restrict ourselves to the case in which M is a linear 
vector space with an inner product ( , ); and we shall write the Poisson 
brackets in the form 

{F,G'} = (VF,J,.VG), 

where Jx is a skew-symmetric linear transformation on M and V -F is the 
gradient of the function F . The gradient is characterized as follows. Differ- 
entiating F{x{t)) along a curve x{t) on M, we have 

If is non-singular, then the Poisson brackets are non-degenerate and 
have locally a canonical system of coordinates. In many problems of physical 
interest, however, the Poisson brackets are degenerate, i.e. det Jx = 0. For 
example, in the study of rigid motions about a fixed point in M^, the Poisson 
bracket is 

= (VF, XX VG). (3.1) 

The operator Jx is defined by Jyy = x x v; hence ker( Jx) = Mx. 

The bracket (3.1) vanishes for all regular functions G whenever F is spher- 
ically symmetric. Such a function F is called a Casimir. It is invariant under 
any Hamiltonian flow generated by these brackets. 

Any Poisson bracket on an odd-dimensional manifold must be degenerate 
and therefore have Casimirs. The bracket (3.1) is an example of a non- 
canonical Poisson bracket. 

The formalism of Poisson brackets and Hamiltonian flows can be extended 
to inflnite dimensions, for example, in the study of continuum mechanics, 
though a number of technical difficulties arise. In particular, Poisson struc- 
tures play a useful role in the theory of the Euler equations for an incompress- 
ible fluid. Two important such brackets are the Poisson bracket introduced 
by Arnold [3, 4, 5] in his study of incompressible fluids on flxed domains, and 
the Poisson bracket implicit in Zakharov's fundamental discovery [18] of the 
Hamiltonian structure of the Euler equations of gravity waves. 

3.1 Arnold's Poisson brackets 

Arnold observed that Euler's equations for an incompressible fluid in a flxed 
domain T) are directly analogous to his equations for rigid body motion, and 
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that they have a Hamiltonian structure with the Hamiltonian and Poisson 
brackets given respectively by 

H = jjjU-wd?^, (3.2) 

V 

and 

{F,G}= [[[ fcurl vx (3.3) 



5v V 6v 



V 



Here, F and G are functionals on with gradients in €„. The gradient 
of F is SF/Sv, the Euler-Lagrange derivative of F with respect to v. For 
example, ff- = v. The operator in this case is 

JvW = Po-(curl V X w). 

Let us show that (2.1) are the Hamiltonian equations generated by (3.2) 
and (3.3). We have 

{H,F} = jjj^- (curl V X ^)rf3x = yjy^ . (v X curl v) d=^x. 

V V 

The Hamiltonian flow F = {H,F} implies that 

6F 

— ■ (vt + (curl v) X v) d^x = (3.4) 
ov 

for all admissible F on 

All linear functionals of the form F^{\) = (w, v) are admissible, and 
the gradient of F^ is the vector w. Therefore + (curl v) x v belongs to 
C-^ = Ct,. Hence it is a gradient, and 

Vt + (curl v) X V = V / 

for some function /. The Euler momentum equations (2.1) follow from this 
and the vector identity 

(curl v) X V = (v ■ V)v - i V |vp, (3.5) 
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if we let / = — p + I V |vp. 

Just as in the case of rigid motion, the Arnold bracket is degenerate. This 
degeneracy is related to the action of the (formal) group of volume preserving 
diffeomorphisms acting on P. Arnold's Poisson bracket is an example of a 
Lie-Poisson bracket. 



3.2 Zakharov's Poisson brackets 



In 1968, Zakharov made a striking observation: Euler's equations for irrota- 
tional gravity waves have a canonical Hamiltonian structure. The Hamilto- 
nian (in non-dimensional variables) is 



H =- 
2 




(V^)2rf3x+-A // C'(x,y,t)rfV 



The Poisson brackets implicit in Zakharov's observation are the canonical 
brackets 

f5F5G 5FSG\ 



the Hamiltonian flow is then the canonical flow 



Ct 



6ip ' 



5H_ 



The Hamiltonian H is regarded as a functional of (^, Q where C, = 
({x,y,t) is the height of the free surface, and ip = ip\s is the trace of the 
harmonic function ip on the free surface, with ^9,^ = on the bottom. The 
evolution takes place in the space of harmonic functions on V. 

Zakharov's result is verified by calculating the gradients of H with respect 
to C and ip. Now 



d 

de 



e=0 



5Cd\ 



where V denotes V (p I 



By identification, 
6H 1 
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Similarly, 



d 

de 



e=0 




V 




dip 
5ip—-dS 



dip 
Sip-dS, 



V 



since p is harmonic in D and p^ = Q on the bottom. 
On the free surface 



p)ydS = Vp 



1 + Q + 



so 



5H 

5p 



P^z - PxCx - PyCy 



The free boundary equations (2.2) and (2.3) are thus precisely the Hamilto- 
nian equations for this system. 

Remark. The effects of surface tension can be obtained by simply adding 
the boundary integral 

a ff dS 



to the Hamiltonian, where a is the coefficient of surface tension, and dS is 
the element of surface area on the free surface S. The inclusion of surface 
tension leads to an additional term in the Bernoulli equation; when the free 
surface is a graph z = ({x,y,t), it is 



ipt + ^1 Vv?P + g z = a div . '^^ =. 



vc = (Cx.,g. 



The potential energy can also be written as the integral of the gravitational 
potential over the fluid domain, so that the Hamiltonian for gravity waves 
including the effects of surface tension is 



H 




V 



{Vpf 



d^x. + a 



dS, 



(3.6) 
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where t/+(x) is the gravitational potential, truncated in such a way that 
the integral over the unbounded domain T) converges. When the fluid is a 
horizontal layer and the gravity field is constant in the negative z direction, 
we take U+ = {z — 1)+, where 2+ denotes the function given by z when 
2; > and by when 2; < 0. The factor g has been absorbed into the pure 
parameter A. 



4 Free boundary flows with vorticity. 

Free boundary value fiows with vorticity, with both gravitational forces and 
surface tension included, are generated by the Hamiltonian 

H = Jjj ££x + a jj dS, £ = ^ + Af/+(x) (4.1) 

V E 

The corresponding Poisson brackets are [15] 



V 



where S is the free boundary and dS is the element of surface area on E. 

Admissible functionals are regarded as functions of v and E, the free 
boundary of T>, and their gradients are defined implicitly by the relation 



, , , — • 5v rf^x + / / — 5S dS. 

e=0 I I I OV I I 01^ 




V 



Variations with respect to the free surface are restricted to normal variations, 
in a sense explained below. Admissible functionals F are those for which 
SF/6\ is a divergence free vector field. We require that jj^5T,dS = 0, 
refiecting the fact that only volume preserving variations are allowed. This 
means that the gradient of a functional with respect to S is determined only 
up to a constant. 
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Let C.(i{T>) be the space of divergence free vector fields on V. Let Pi 
and P2 be the orthogonal projections defined by Piv = w and P2V = V in 
the Weyl-Hodge decomposition. 

Lemma 4.1. Let 

F(v, S) = jjj J^(v, x)rf^x + a jjdS 
V s 

6e an admissible functional. Then 

6F „ 5F 



ow ov 



The gradients with respect to ip and S lie in >C^(S, dS) and are given by 

mod constant. 



SF _ 6F 

where k is the mean curvature function on S. 



■V, — = T[y, x) + OK 



Proof. Applying the Weyl-Hodge decomposition to both 6v and 5F/5v we 
obtain 



/ J J J J J J 

By the uniqueness of the Weyl-Hodge decomposition, we may conclude 

5F _ p5F 5F _ p5F 

Since 5F/5\ is divergence free, we have, by the divergence theorem, 

P2^-5Vv^c/=^x = JJJ ^■\/6ipd^x = JJ^.uSipdS, 

V V dV 

and the second relation follows. 

Let Eg be a one parameter family of surfaces parameterized by a vector 
valued map 

X(ti, u, e) = Xo(u, v) +e5'L N(u, t;). 
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where N is the normal vector field to S. For e sufficiently small, the sym- 
metric difference T>^AV of the domains bounded respectively by Eg and S is 
contained in a tubular neighborhood of S. In this neighborhood, the volume 
element of the 3-space can be written as d'^x = drdS where dS is the area 
element on S and dr corresponds to the normal coordinate in the tubular 
neighborhood. We get 



Sill J^(v, x) d^x = lim i 





J^(v, x) d^x 



J j (^lim ^ ^ J^(v, x) rfr^ = JJ ^)^^ 



s s 
On the other hand, by classical differential geometry. 



S JJ dS = JJ t^^^ dS, 
s s 

where k is the mean curvature function on S. This completes the proof of 
Lemma 4.1. □ 



Let us derive the equations of motion from the Hamiltonian structure. 
We have 

6H 6H 6H ^ 

— — = V, — — = V ■ I/, — = C + CTK. 

ov dip oS E 



From F = {H, F}, we get 




SF_ 
5^ 



vt d^x. + 



S+ dS 



V 




y X (curl v)) ■ ^ (i^x 



V 

Since 



5F SF^ 
v-i.^-(£: + aK)^) dS. (4.3) 
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and 5F/5\ is divergence free, we get from (4.3), using functionals for which 

5F/5ip = 0, 



Vt + (curl v) X V = V (— p + S), = V ■ u 



s 



The boundary condition on the bottom is v ■ = 0, where u is the outward 
normaL 

The first equation, together with (3.5) imply that Ap = — div (v • V)v. 
Substituting the two equations above into (4.3), we obtain 

W p-^d^x- [[ aJ-f- dS = 



V 



5v J J Sif 



for all admissible functionals F. Applying the divergence theorem to the 
integral over V we obtain 

6F 

(p — an)^— dS = 0, 

dip 
s 

for all admissible functionals F. But 



and therefore 



p 



= an + constant. (4-4) 

s 

Thus the Hamiltonian approach yields the dynamic conditions on the free 
boundary in the case of surface tension. [7, 15, 12]. 

Remark. In the general theory one considers normal variations of the free 
surface, whereas in the theory of gravity waves on a free surface over a hori- 
zontal bottom, it is customary to use the height of the free surface, (. More 
generally, if the surface S is a graph over a fixed manifold A4, we may rep- 
resent S by a "height" function ( defined on M. In that case we refer to 6( 
as the "vertical" variation and 5S as the "normal" variation. 

Proposition 4.2. Let (5S and denote the normal and vertical variations 
of a surface S in the case when T, is a graph over a fixed manifold. Let E be 
given in local coordinates by (p = 0, where (j) = z — (. Then S( = | V<^| 5E. 
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Proof. Let X : f/ i— > be a local embedding of S in M^; and let be a 
one parameter family of embeddings, with Xq = X. Then 



de 



e=0 



Let S be defined by 
d 



5T. 



de 



0, (f> = z — (. Then X^ = {u, v, (ei^u, v); and 

V0 S( 













V 






e=0 









□ 



5 Variational principles for traveling waves 

The Hamiltonian structure of the equations for gravity waves can be used to 
obtain variational principles for traveling waves - waves of constant speed and 
shape. Such a wave is a stationary solution of the Hamiltonian system in a 
Galilean frame moving with the wave; thus the wave is a critical point for the 
Hamiltonian, computed in such a reference frame. We apply the method here 
to the general case of gravity waves on a horizontal surface. The variational 
principle for irrotational flows given below appears to be new. 

A variational approach, if successful, would permit a global treatment of 
the existence of traveling waves by the direct methods of the calculus of vari- 
ation; but so far, the existence of traveling waves for potential flows of low 
amplitude have been proved by perturbation methods. The first existence 
theorems were given independently for periodic wave trains by Levi-Civita 
[14] and Struik [17] in the case of finite depth. The existence of the solitary 
wave, which is a more difficult problem, was first proved by Friedrichs and 
Hyers [10], since the bifurcation problem in this case is a singular pertur- 
bation problem (see the discussion by Sattinger [16]). These authors used 
conformal mapping techniques. A dynamical systems approach to the exis- 
tence of traveling waves has been developed by Kirchgassner [13]; Amick and 
Toland [2] have shown that periodic wave trains tend to a solitary wave in 
the limit as the period tends to infinity. 

In the direct method, one first uses compactness properties of the func- 
tional to obtain a minimum from a minimizing sequence. In general, this 
guarantees only a weak solution of the associated Euler-Lagrange equations. 



15 



In many cases, these are elliptic equations, and it is possible to prove sufficient 
regularity of the weak solution to show that in fact it is a classical solution to 
the problem. (See Alt and Caffarelli [1] for functionals of Friedrichs' type.) 
For the present, we simply indicate the method for the problems discussed 
here in the theorems below. 

Theorem 5.1. Euler's equations for gravity waves are the Euler- Lagrange 
equations for the functional 



rf^x, (5.1) 



where 



0, y<0; 

y, y>o. 



T^C = {{x, y) : -oo <x<oo, ^ <y <l + Ci^)] ; 
and the minimum is taken over all functions ip for which 

[{^^ - if + (^J] rf^X < +00. 



V 



If{ip,Q is a local minimum ofTi, then is harmonic on the interior ofD,^; 
if C, is and ip G H'^{T>), then the kinematic and Bernoulli equations hold 
on the free surface. 

Remark. The Hamiltonian (5.1) is the renormalization of the Hamiltonian 
in the moving frame. By carrying out the integration in y we obtain 



jj{y-l)^d''^=\j\'dx- 



V 

thus Ti can also be written 



oo 

2 



C (x) dx. 



V 



If (/? G H^{T>) and ^ G C^, then V has an LP' trace on the boundary y = (, 
and Stokes theorem applies. 
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Proof. Let (</?, () be a minimizer of Ti and suppose that ( is and (f G 
H'^{T>). Let (^g) be a one parameter family of admissible functions and 
denote the corresponding domains by P^. By the calculations in §3.2 we have 



de e=o 
Aip5ip (fx + y (pu5^ ds 

V dV 



+00 



0, 



5C, dx 



for all admissible 6ip, 6(. 

Since the bottom is fixed, ip^ = on y = 0. We first restrict ourselves to 
variations for which = 5(p\^ = 0. Then the double integral must vanish for 
a set of variations 6(p which are dense in it follows that (f is harmonic 

in the interior of V. As before, ip^ds = V f ■ {—(x, l)dx = {ipy — (x^x)dx; 
and so 



6n 



+00 



1 



AC MC + (v^j/ - fxCx) Sip 



dx. 



Setting first 6( = and letting S(p vary on S, we obtain the kinematic 
equation on the free surface. Therefore the second term always vanishes. 
Now allowing S( to vary, we see that Bernoulli's equation holds on S. □ 



6 A variational problem with constraint 

Whereas Friedrich's paper shows that Bernoulli's equation is not obtained 
when the functional J is minimized with respect to the stream function, 
Constantin et. al. showed in [7] that traveling gravity waves in the rota- 
tional case are obtained as extremals of a variational problem for the stream 
function with constraints. The existence of traveling water waves with vor- 
ticity was established in [8] for the periodic case. In a recent PhD thesis at 
Brown University, V. Hur [11] has constructed solitary waves with non-zero 
vorticity. Some of their qualitative properties were investigated in [6]. 
In the irrotational case we have 
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Theorem 6.1. Define the set of admissible functions K, = {ip,C} with the 
following properties 



+00 p+00 

2 



i) / ({x)dx = m; / C, dx < 00; 

J —00 J —00 

ii) V(a;,0) = 0; ^/'(x, 1 + C(x)) = 1, 
Hi) jj [^1 + {ipy - 1)^] rf^x < +00. 

V 

Consider the variational problem 

//^ [(VV^)^ - 1] ci^x 



A = inf ■ 



J —00 ^ 



Let [ip, Q be a minimizer in K, of the above variational principle. Then ip is 
harmonic in the interior ofD. If C, is C^, andijj G H'^{T>), then the Bernoulli 
equation is satisfied on the free surface if) = 1. Hence minima of the above 
variational problem provide an irrotational flow for the gravity wave problem. 

Proof. Let {ip^Q be a minimizer, and let ipsiCe be a family of admissible 
functions with ipQ = ijj and Co = C- Then J{e) > and J(0) = 0, where 

r r r+co 

J{e) = jj [(V^,)2 - 1] rf^x - A y C dx. 
Then 6J{6ip, 6() = for all admissible variations, where 

r r r+00 

SJ= 2 W ■ V 6tljd^x+ / [{Vtlj)^ - 1 - 2AC] SCdx 

= -2 JJ Atp6ijd'^^ + 2 ^ 6ij^lj^ds + J [(V^^)^ - 1 - 2AC] ^C^^a;. 

The integral over the bottom of the flow domain vanishes, since = 
there. On the free surface (see S{2), p. 65 [9]) 

(5^ + ijy6C = 0. 
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This follows immediately by differentiating the relation ips{x, 1 + Cei^)) = 1 
with respect to e and setting e equal to zero. Similarly, differentiating the 
expression ip{x, 1 + C{^)) = 1 with respect to x we find that ip^/ipy = —(x] 
hence 



Hence 5 J reduces to 



r r r+oo 

6J = -2 JJ Atlj6tljd'^:>c- J [{V tpf + 1 + 2XC] 6C dx. (6.1) 



First restrict the variations to fixed domains, 6( = 0, and the first integral 
must vanish for all variations 6ip which vanish on dV. Hence ijj is harmonic 
in the interior of P, and the double integral vanishes. 

We next consider variations of the domain. Since J (^dx = m, for all 
variations, we have J 6( dx = 0; then the condition 



/ + CO 
{{V + 1 + 2XC)6C dx = 
-oo 



for all such S( implies that the integrand is a constant. We therefore have 
(V + 2X( + 1 = C = const, on the line; letting x —>■ oo and noting that 
C — »■ while (V ^ 1 we see that C = 2, and the Bernoulli equation is 
satisfied. □ 
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